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Graphene is a truly two-dimensional material with ex- 
ceptional electronic, mechanical, and optical properties. 
As such, it consists of surface only and can be probed by 
the well developed surface-science techniques as, e.g., 
scanning tunneling microscopy. This method bridges the 
gap between the surface science community and the elec- 
tronic device community and might lead to novel com- 
bined approaches. 

Here, I review some of the scanning tunneling mi- 
croscopy (STM) and spectroscopy (STS) experiments 
on monolayer graphene samples. I will concentrate on 
graphene samples deposited on insulating substrates, 
since these are related to graphene device concepts. In 
particular, I will discuss the morphology of graphene on 
SiC>2 and other emerging substrates, some nanomechan- 
ical manipulation experiments using STM, and spectro- 
scopic results. The latter can map the disorder potentials 
as well as the interaction of the electrons with the dis- 
order which is most pronounced in the quantum Hall 
regime. 




3D-representation of a scanning tunneling microscopy image 
of graphene; the atomic resolution is displayed as color code, 
while the rippling is shown in 3D. (courtesy of M. Pratzer, 
RWTH Aachen). 
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1 Introduction Since its discovery in 2004 [1|, 
Graphene has become one of the most investigated ma- 
terials worldwide with more than 3500 publications per 
year. This is partly due to the exceptional band structure 
of graphene which might mimic Dirac physics in table-top 
experiments EJO, but increasingly due to the applicational 
prospects with respect to, e.g., transparent electrodes [4|, 
high frequency transistors (3), optoelectronic devices [6| 
and sequencers for individual DNA strands (7). Further 
applications in spintronics OH or quantum computation 



1 10 1 are envisioned. 

Particular advantages of graphene are the abundance of 
its constituent atoms, its inertness, its very high break- 
ing strength, its high room-temperture mobility, and its 
simple preparation. Moreover, as a 2D material it con- 
sists of surface only and can be characterized in depth 
by the well developed surface science techniques such 
as angular resolved photoelectron spectroscopy (ARPES) 
ifTTI or scanning tunneling microscopy (STM). The latter 
technique reveals precise information on morphology and 
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electronic structure. In short, graphene offers the unique 
possibility to combine precise knowledge of a device down 
to the atomic scale with functionality. This feature article 
gives a short summary of the major achievements obtained 
with STM so far. It concentrates on graphene on insulators, 
since metal substrates are limited with respect to electronic 
applications due to the shortcircuit created by the metal. 
As usual, only a selection of results can be presented. 



2 Scanning tunneling microscopy In scanning 
tunneling microscopy, a sharp metallic tip, mostly end- 
ing in a single atom, is positioned 3 — 8 A above a con- 
ducting surface. The surface is usually prepared in ultra 
high vacuum (UHV) in order to be atomically smooth. 
The high inertness of graphene allows preparation out- 
side UHV. Subsequent annealing to, e.g., 400° C, is suf- 
ficient to get rid of all unwanted adsorbates leaving an 
atomically smooth graphene surface as presented, e.g., 
in the introductory figure. A voltage V is applied be- 
tween the tip and the conducting surface and the result- 
ing tunneling current / is measured. / depends exponen- 
tially on the distance between surface and tip Az accord- 
ing to I(z) oc e~ aAz . A good approximation for a is 
a = yJ(Am ■ (<? s + <f> t - e\V\)/h) with Planck's constant 
h, electron mass m, electron charge e and work functions 
of tip and sample <P t and <P S , respectively. A good estimate 
is a ~ 2.1/A. 

The tip is positioned with respect to the sample using 
piezoelectric elements. All three directions x, y, and z 
can be changed with sub-pm precision [ 14- 1 . For STM, the 
tip is scanned in x and y direction and the tunneling cur- 
rent is kept constant by a feedback mechanism adjusting 
z. The resulting z(x, y) is plotted and called constant- 
current image. It represents, to first order, a contour of 
constant integrated local density of states of the substrate 
LDOS(x,y, z,E), where the integration has to be taken 
between the Fermi levels of sample and tip to be adjusted 
by V [12]. The central position of the very last atom of the 
tip is given by (x, y, z) and E is the energy. Such images 
are often called topography of the sample indicating that 
corrugations of the atomic positions dominate the contour. 
Differentiating / with respect to V (at low V with respect 
to $ B (t)) gives direct access to the LDOS according to: 

dI/dV{x, y, z, V) oc LDOS{x, y, z, E) = (1) 

^2\^(x,y,z)\ 2 -S(E-E). 

Thereby, ^ are the single-particle wave functions of 
the substrate at energy E and E = eV. Of course, 
this requires that the system is adequately described by 
independent single-particle wave functions. Moreover, 
an s-type symmetry of the orbitals of the last atom is 
necessary, in principle, but Chen has shown that the 
model remains largely correct even if higher orbital mo- 
menta are contributing [EQ. In real experiments, the 6- 
function has to be replaced by an energy resolution func- 



tion with approximate full width of half-maximum of 
SE « ^(3.3 ■ k B T) 2 + (1.8 • eV mod ) 2 . T is the temper- 
ature (fc: Boltzmann's constant) and V mo a is a modulation 
amplitude used to detect dl/dV by lock-in technique. The 
resulting dI/dV(x, y) recorded at constant V and z(x, y)- 
values resulting from constant / is called LDOS'-image. 
Thus, STM can measure atomic structure with sub-pm 
resolution and electronic structure (LDOS) with sub-meV 
resolution [14|. The energy resolution makes STM com- 
plementary to the transmission electron microscope (TEM) 
which reveals atomically well defined structural informa- 
tion partly with chemical specifity [15], but not the LDOS 
down to the meV scale. 

2.1 Graphene Graphene is a two-dimensional (2D) 
honeycomb lattice of carbon atoms, which are bound in 
sp2+Pz configuration leading to a -bonds and 7r-bonds, re- 
spectively [16]. As a 2D material, it should be unstable at 
T > K 07], but it has been found theoretically that a 
coupling of out-of-plane and in-plane phonons could sta- 
bilize graphene without a support leading to a rippled mor- 
phology ifTSI . Such rippling with wave length of about 10 
nm and amplitude of about 1 nm has indeed been found 
experimentally, if the graphene is only partly supported at 
the edges, e.g. around a hole of a TEM grid, and freely 
suspended in between 1191 . It is still debated, if the experi- 
mentally found rippling is fundamental or if it is caused by 
preparation history ll20l . 

The electronic structure of graphene at Ep is governed 
by the bonding and antibonding parts of the 7r-band. These 
bands touch at the corners of the hexagonal Brillouin zone 
(BZ) due to the fact that the unit cell exhibits two atoms, 
which are attributed to two different sublattices. The en- 
ergy dispersion is, to very good approximation, lin- 
ear around Ep leading to conical E(k) planes: E(k) = 
ftvpj\k\ (fc: wave vector, vp> ~ 10 6 m/s: Dirac velocity). In 
first order, the bands are also electron-hole symmetric, i.e. 
a cone is opening upwards into the unoccupied states and 
downwards into the occupied states from the six points at 
the corners of the BZ. The energy at the crossing points is 
called Dirac energy Ep, since it can be moved away from 
Ep by charging the graphene layer, either using dopants or 
a gate voltage. 

In addition, the electronic wave functions show a sublattice 
texture, which changes continuously, if one moves around 
the cone 0. Since the sublattice textures at opposite sides 
of the cone are orthogonal, backscattering of electrons is 
suppressed OTTl . This is regarded as a major reason for the 
high electron mobility of up to /i ~ 10 m 2 /Vs found close 
to room temperature [22] (for comparison: mobility of Si at 
300 K: /j, ~ 0.1 m 2 /Vs). The linear band structure and the 
sublattice texture renders the electrons in graphene equiv- 
alent to chiral massless particles at the velocity of light al- 
beit with the reduced velocity vp>. Thus, graphene raises 
hope that relativistic effects as Zitterbewegung, Klein tun- 
neling or supercritical charge might be observed in table- 
top experiments |[3ll23l. 
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Figure 1 (a) STM image of graphene deposited on Si0 2 by the scotch tape method, V = 1 V, I = 0.2 nA, T = 300 K; 
(b) optical microscope image of a monolayer graphene sample contacted by In microsoldering [28 50 1 ; (c) same as (a), 
but on another sample, V = 0.4 V, I = 0.2 nA, T = 300 K; (d) scanning electron microscopy image of the W tip used 
for the AFM measurements shown in (f), nominal curvature radius at the apex is indicated; (e) oscillation amplitude of the 
cantilever as a function of tip-sample-distance, the tip-sample distance used for recording (f) is marked by an arrow; (f) 
Tapping-mode AFM-image of the SiC>2 substrate used for the preparation of (a) and (b), Af — —100 Hz, force constant: 
47 N/m, oscillation amplitude: A = 18 nm [36 1; (g) same as (a) and (c), but measured by another group, V = —0.3 V, 
/ = 0.04 nA, T = 300 K; (h) non-contact AFM image of the Si0 2 used for the preparation of (f), Af = -20 Hz, A = 5 
nm [37]; all images execept (f) are recorded in ultrahigh vacuum, (courtesy of W. G. Cullen, University of Maryland for 
(g),(h)) 



In magnetic B-field, the linear dispersion of graphene is di- 
vided into Landau and spin levels. The energy of the Lan- 
dau levels is given by E N = ±v D ■ y / 2eh ■ \N\B (N £ Z). 
The square-root dependence on \B\ leads to a relatively 
large separation of the central Landau level Eq, which, 
protected by symmetry, is located at Ed, and the two sur- 
rounding Landau levels. An energy distance of AE = 0.1 
eV can be reached already at B = 5 T, which is partly re- 



sponsible for the fact that quantum Hall plateaus can be ob- 
served at 300 K [24|. Consequently, it is attractive to study 
quantum Hall physics, which is driven by interaction of 
electrons with disorder and, thus, by inhomogeneous elec- 
tron distribution [25), using a local scale method like STM 
on graphene. 

Also the more intricate fractional quantum Hall effect [25ll 
is significantly more stable in graphene than in the conven- 
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tionally used III-V semiconductors [26,27]. This can be 
rationalized, since the energy gaps, which cause the addi- 
tional plateaus in the Hall curve at particular ratios between 
electron density n and density of magnetic flux units till, 
are driven by electron-electron interactions [25 1. These in- 
teractions are much more intense within a single atomic 
layer surrounded by a dielectric material with e = 1 than in 
a 2D electron system in III-V semiconductors with thick- 
ness 10 nm and surrounded by e = 13. The energy gaps 
E gap at v := n/riLL = 1/3 of freely suspended graphene 
were indeed found to be -E gap = 4 meV at B = 20 T fl26l 
and even after depositing the graphene on BN(0001), they 
were E sap = 2 meV at B = 12 T 11271 . Thus, studying 
fractional quantum Hall physics on the nm scale might be 
within reach on graphene. 

Notice that each Landau level of graphene is fourfold de- 
generate due to usual spin degeneracy and due to the two 
non-equivalent corners of the Brillouin zone K and K', 
which leads to so-called valley degeneracy. 

3 Graphene Morphology The first graphene sam- 
ples of one monolayer height have been prepared by 
the scotch-tape method on SiC>2 |lj, but shortly after re- 
searchers also managed to show single-layer behavior for 
the sublimation method applied to SiC(OOOl) 11291 - Since 
the later samples cover a complete area, they are easier to 
find by STM and, thus, have been imaged rather early in 
ultra high vacuum [ 30 1 . It has been found that single layer 
graphene exhibits a y3 x \/3 reconstruction 1 30 1 being 
located on a strongly bound carbon layer in sp3 configura- 
tion inn . 

But also graphene flakes deposited on Si02, which exhibit 
higher mobility and mostly better defined quantum prop- 
erties [2 , 32 33 1 have been imaged by STM. They exhibit a 
rippling on the 10 nm length scale with amplitudes of about 
1 nm Q4II35I . Figure [TJ a), (c), and (g) show three STM 
images of typical graphene flakes prepared by the scotch 
tape method in comparison with atomic force microscopy 
(AFM) measurements of the respective Si02 substrates 
as shown in Fig. |TJf) and (h) 03611371 . Partly, the length 
scales between substrate and graphene corrugation coin- 
cide (Fig. [TJc)/(f) and (g)/(h)), but partly the length scale 
of the graphene rippling is a factor of 3-4 smaller than the 
length scale of corrugation on Si02 (Fig. [Tfa)/(f)). Since 
the latter might be an artifact of reduced lateral resolution 
of the AFM measurement, the authors used a particularly 
sharp tip as shown in Fig. [TJd) and controlled the con- 
tact depth by oscillation amplitude vs. distance curves as 
shown in Fig. [JJe). This provides a lateral resolution of 1 
nm much shorter than the observed length scale of rippling. 
A detailed analysis using Fourier transformation and auto- 
correlation functions of the real space images revealed that 
the latter type of samples, which are found in 80 % of the 
preparations reported in l36ll . are partly freely suspended 
above the substrate and exhibit a slightly preferential wave 
length of 15 nm, which surprisingly is very similar to the 



rippling found on freely suspended graphene flakes by 
TEM (19] . However, the image within the abstract, which 
is taken on the same sample, shows that rippling on all 
length scales down to about 2 nm is present. (36]. 
The other samples (Fig. [TJa)/(f) and (g)/(h)) show a nearly 
perfect conformation of the graphene flake to the substrate 
. This is to be expected due to the attractive van-der-Waals 
forces between the substrate and the graphene outperform- 
ing the elastic restoring forces of the graphene [37]. The 
latter are fully characterized by the elastic modulus mea- 
sured to be E = 340 N/m |38l . 

Obviously, the details of the morphology of graphene on 
Si02 depend on the uncontrolled details of the prepara- 
tion procedure, which are also apparent, e.g., by mobility 
values varying by a factor of 100 l2l [39l . So far, it is un- 
clear, in how far the rippling is responsible for the disorder 
potential and the mobility of graphene [40|. However, de- 
positing graphene on an intermediate BN layer decreases 
the rippling down to a value of about 0.01 nm PT1 and, at 
the same time, reduces the disorder potential by more than 
a factor of 10 leading to strongly improved mobility ll26l . 



4 Nanodrums The partly freely suspended graphene 
samples on Si02 can be used for nanomechanical manipu- 
lation of the one atomic layer thick membrane. Therefore, 
the tip of the STM is used fl4l . whose forces on a local 
graphene area can overcome the van-der-Waals forces of 
the substrate. Figure |2ja)-(c) show that valleys of the rip- 
pled graphene can indeed be lifted by the tip forces. The 
central valley is transformed into a hill by decreasing the 
distance between tip and substrate and, thereby, increas- 
ing the respective van-der-Waals force. Figure |2|d) depicts 
that the corresponding images are dynamic, since the force 
between valley and tip changes during the imaging proce- 
dure. However, the I(z) curves taken within the center of 
a valley can be used to classify the valleys. 10 % of the 
valleys exhibit a tip-sample distance region of about 0.15 
nm, which exhibits a stronger exponential increase of tun- 
neling current than expected from the work function (Fig. 
|2je)). This indicates continuous lifting of the graphene val- 
ley. 50 % of the valleys do not move at all probably because 
they touch the substrate. The remaining 40 % exhibit a hys- 
teretic I(z) curve with jumps in tunneling current by about 
three orders of magnitude as shown in Fig.|2jf). This indi- 
cates bistable behaviour of the membrane. Indeed, adding 
up the elastic restoring potential of graphene and the two 
van-der-Waals potentials of tip and sample reveals a dou- 
ble well potential lfl4ll as shown in Fig. |2jg). This double 
well potential can be additionally tilted by the dielectric 
force (image charge) of the biased tip. Thus, a flipping of 
the membrane between the two potential valleys is possible 
either by applying a tip voltage as shown in Fig. |2|h) and 
(i), or by changing the tip-sample distance and, thereby, 
changing the van-der-Waals force of the tip. 
An interesting aspect of the movable nanomembranes 
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Figure 2 (a) STM image of graphene deposited on SiC>2, V = 1 V, I = 0.1 nA, T = 5 K, atomic resolution is shown 
as color code, while long-range rippling is displayed in 3D; (b) same as (a), but at V = 0.7 V, I = 1 nA; (c) same 
as (a), but at V — 0.4 V, I = 2 nA; small images indicate decreasing tip-surface distance; (d) sketch of the dynamic 
movement of the valley during scanning of the tip (black lines) and the resulting, measured line (red); the directions of 
the four forces acting on the membrane are color-coded; (e) I(z) curves measured on a valley, which moves continuously 
towards and away from the tip; voltage V and stabilization currents I s tab are indicated; the red dashed lines mark the 
decay ln(///o) = 2^2m<j)/Ti (m: electron mass, <j) — 4.85 eV: work function averaged between tip and surface, h: 
Planck's constant, Iq: current at saturation) expected from pure tunneling; (f) I(z) curve measured on a valley which 
moves fry steretic ally; V = 1 V; directions of tip movement are indicated by arrows; (g) diagram showing the different 
potentials acting on the membrane as well as the sum potential ^ S um; ^mcm- elastic restoring potential of the graphene 
valley, ^ i: potential induced by the image charge of the biased tip, $ v dw(tg) : van-der-Waals potential between tip and 
graphene, ^ v dw(sg) : van-der-Waals potential between SiC>2 substrate and graphene; z sg : distance between graphene and 
substrate; (h), (i) sum potentials at the tip voltages indicated and at a tip-substrate distance of 1.46 nm; the arrows mark 
the movement of the membrane induced by switching the voltage between these two values. [ 14 1 



is the presumably extremely large resonance frequency, 
which can be estimated to be 0.5 Thz 11411421 . It is given 
by the small number of contributing atoms (800), its light 
mass (12 u) and the large elastic strength of graphene with 
an elastic modulus of i?2D=340 N/m. The corresponding 
vibrational energy of a movable valley is 2 meV, leading to 
99.9 % ground state occupancy already at T — 4 K and a 
complete suppression of excitations at T = 0.3 K, a tem- 
perature which is meanwhile achieved routinely in many 
STM laboratories |43l. 

In order to estimate the required electric field amplitude 
for exciting at resonance, we measured the oscillating me- 
chanical response of the membrane to an oscillating volt- 
age applied to the tip. Figure |3jb) and (c) show the applied 
voltage and the resulting current response on a non-moving 
area in (b) and on a continuously moving valley in (c). The 
former exhibits only a 90° phase-shifted sinusoidal current 



given by the capacitive cross-talk between tip and graphene 
sample. However, the latter clearly shows an additional, 
non-sinusoidal in-phase signal indicating the movement of 
the sample. The amplitude of this signal can be recalcu- 
lated to a height amplitude of the oscillating membrane 
by using measured I(z) curves on non-moving areas [14|. 
The result is shown for different excitation amplitudes in 
Fig.[3jd) revealing that an amplitude of 1 V/nm is required 
for an oscillation amplitude of 0. 1 nm at the excitation fre- 
quency of 1000 Hz being far below resonance. 
Taking the quality factor of 10 5 recently obtained on larger 
graphene flakes at much larger amplitude |44|, the required 
voltage amplitude would be reduced to 10 4 V/m, which 
should be taken as an orientational value for future experi- 
ments heading for a resonant coupling to these nanomem- 
branes. 
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Figure 3 (a) STM image of graphene deposited on SiC>2, 
V = -0.4 V, I = 2 nA, T = 5 K; tip with applied ac volt- 
age V(t) and resulting current I(t) is sketched; yellow and 
grey dots mark measurement positions for (b) and (c), re- 
spectively; (b) applied ac voltage V(t) — V mo d ■ sin(oj ■ t) 
(blue) and resulting current response I(t) (red); measured 
at the position marked by a yellow dot in (a), tip was sta- 
bilized at V = -0.4 V, I = 2 nA; (c) same as (b), but at 
the position marked by a grey dot in (a); (d) In-phase I(t) 
amplitude at applied frequency cj = 2 ■ ir ■ 1000 rad/s mea- 
sured at the grey dot in (a) as a function of voltage ampli- 
tude Vi no d; top and left scale show electric field amplitude 
E and resulting oscillation amplitude of the membrane Az 



5 Scanning Tunneling Spectroscopy First spec- 
troscopic results have been published for graphene on 
SiC(OOOl), which exhibit a clear Dirac cone dispersion by 
analyzing the wave length dependence of scattering states 
11451 . but the dl/dV curves remain elusive l46l . 
On graphene flakes, Zhang et al. l47l observed a gap 
around the Fermi level which appears to be independent 
from the back gate voltage moving the Dirac point through 
the Fermi level. This gap of about 130 meV is shown in 
Fig.[4ja). It has been related to the fact that tunneling into 
the graphene is preferentially to electronic states at the 
-T-point, but not to the if -point, where the Dirac cone is 
located. Tunneling into the Dirac cone, thus, requires addi- 
tional momentum, which, e.g., can be delivered by exciting 
a phonon. The lowest energy of a phonon with sufficient 
momentum is the out-of-plane accoustical phonon at the 
K point, which, indeed, has an energy of about 65 meV 
|48| naturally explaining the gap of ±65 meV. Additional 
evidence comes from the stronger decay of the tunneling 
current I with tip-surface distance z within the gap (Fig. 
|4|b)), where direct tunneling into If -point electrons is 
required. A fit by the usual formula to the decay length 

A -1 = 2 • . 2m<P/Ti + kf, (m: electron mass, (f> = 4.85 




sample bias [mV] 



Figure 4 (a) dl/dV curves recorded on the same posi- 
tion of graphene deposited on SiC>2 at different back gate 
voltage Vjgatc, VJjate is given on the right of the curves, the 
arrow marks the Dirac point, Ktab = 0.5 V, J s tab = 0.1 
nA, T = 4.8 K; (b) inverse decay length 1/A deduced from 
fitting I(z) curves by I(z) oc exp (— z/X) as a function of 
sample bias V; I s tab — VI (30 GQ); symbols for measure- 
ments on graphene and Au(lll) are indicated B7l (cour- 
tesy of M. Crommie, UC Berkeley); (c) dl/dV curves 
recorded on the same position of graphene deposited on 
SiC>2 at different / sta b as marked, V sta b — 0.5 V, T = 5 K 

ED. 



eV: work function average tip/surface, h: Planck's con- 
stant, fc| | : wave vector parallel to the surface) indeed nicely 
fits the data. 

However, other authors [49 1 did not observe the gap or do 
observe the gap at low current only. The latter is shown in 
Fig. |4| C X where a gap of 90 meV is observed at stabiliza- 
tion currents below 100 pA, but closes at higher current 
[ 50 1 . Also measurements of graphene on graphite did not 
observe the gap, but a linearly increasing dI/dV(V) away 
from the Dirac point (see Fig.|7ja)) ||5T). Thus, it is likely 
that details of the tip-sample geometry determine, if direct 
tunneling from the tip to A'-point electrons is suppressed. 

5.1 Standing Waves and Potential Fluctuations 

STS can be used to map standing waves of Dirac electrons 
scattered by potential disturbances. Zhang et al. [52 1 have 
found such standing waves on graphene flakes on SiC>2 in 
dl/dV images as shown in Fig.|5jb). The Fourier transfor- 
mation of the real space data (Fig. |5|c)) revealed circular 
structures indicating a dominating wave vector value \k\ 
resulting from intracone scattering. Plotting this value as a 
function of bias voltage revealed a gapped linear 
dispersion in accordance with the gap described in Fig. 
[4] But the steepness of the linear dispersion results in a 
Dirac velocity of vjj ~ 1.4 • 10 6 m/s, which is about 
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Figure 5 (a) STM image of graphene on 285 nm Si0 2 , T = 4.8 K; (b) dl/dV image of the same area, V = -0.75 V, 
/ = 0.08 nA, Vg a t = 60 V; crosses mark the centers of standing waves; (c) Fourier transformation of dl/dV images 
obtained at V = 0.35 V, I = 0.05 nA (top), V = 0.6 V, J = 0.06 nA (middle), and V = 0.85 V, I = 0.07 nA (bottom), 
^gate = 15 V ; (d) dl/dV image of the same area as (a), (b), but recorded close to the Dirac point, V = —0.29 V, 
/ = 0.025 nA, Vg a tc = 15 V; crosses at positions as in (b) [52|; (e) Dirac point map deduced from the minimum in dl/dV 
curves of graphene on 285 nm Si02, T = 4.5 K; (f) same as (e) for graphene on 14 nm hexagonal BN on top of 285 nm 
Si02 ; (g) histogram of the Dirac point values obtained in (e) and (f); (g) histogram of the height values obtained for 
graphene on BN/Si02 and Si02 PTlk ((a)— (d): courtesy of M. Crommie, UC Berkeley, (e)— (h): courtesy of B. LeRoy, 
University of Arizona, Tuscon). 



35 % larger than expected [3]. Since the authors could 
change the gate voltage, they could exclude an influence of 
tip induced band bending [53,54] on this large value. The 
circular Fourier transformation pattern indicates intracone 
backscattering, which requires short range scatterers on 
monolayer graphene ETI . A mapping of the lateral Dirac 
point distribution, performed as a dl/dV image slightly be- 
low the Dirac point (Fig. [5jd)), revealed that the centers of 
the standing wave patterns induce a lowering of the Dirac 
point, which is interpreted as evidence of a charged donor 
character of the scatterers. 

Standing wave patterns on bilayer material on SiC(0001) 
also exhibit contributions from intracone as well as from 
intercone scattering processes l45ll with the Fourier trans- 
formations of the intracone part revealing a linear dis- 
persion gving a Dirac velocity of 1.0 • 10 6 m/s. Fourier 
transformations of standing wave patterns from monolayer 
graphene on SiC(0001) were interpreted as evidence for 
the absence of intracone backscattering |55"1 . 
The Dirac cone mapping in Fig. (5jd) revealed potential 
fluctuations by about ± 30 meV on a 10—20 nm length 



scale. Xue et al. ATI observed a much larger fluctuation 
of ±150 meV on length scales of 3—10 nm as shown 
in Fig. |5je) again pointing to the difference of nominally 
identically prepared samples. A strong reduction of poten- 
tial fluctuations is achieved by putting an intermediate BN 
layer between the graphene and the Si02- This reduces 
the potential fluctuations to ±5 meV at a length scale of 
50—100 nm as shown in Fig. |5jf) [41 1. The difference be- 
tween the substrates is emphasized by comparing the his- 
tograms of potential values displayed in Fig. |5jg). Inter- 
estingly, the difference of histograms is very similar to the 
difference of histograms of the topographical corrugation 
of graphene on the two substrates. This is evident by com- 
paring Fig. [SJg) and (h). Notice that other authors found 
larger potential corrugations of ±15 meV on graphene on 
BN/Si0 2 & 

5.2 Edge States A very intriguing property of 
graphene is that unreconstructed zig-zag edges, which are 
singly terminated, should exhibit a magnetic edge state, 
which, moreover, is antiferromagnetically coupled to the 
zig-zag edges terminated by the other sublattice (57]. In- 
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Figure 6 (a) 3D representation of an STM image of a 
graphene nanoribbon of width 19.5 nm having an (8,1) 
edge, V = 0.3 V, / = 0.06 nA, T — 7 K, the arrow 
with dots marks the positions where the spectra in (b) are 
recorded; (b) dl/dV curves obtained on the nanoribbon 
shown in (a) at the dots marked; distance from the edge is 
indicated in A and peaks are highlighted by dashed lines; 
curves are offset for clarity; inset shows the spectrum at the 
edge of a smaller ribbon with the energy distance of peaks 
A marked, V a t a h = 0.15 V, istab = 0.05 nA; (c) peak 
distance A as a function of nanoribbon width (symbols) in 
comparison with predictions from a Hubbard model calcu- 
lation (pink area); the angles are measured with respect to 
the zig-zag direction ll6D (courtesy of M. Crommie, UC 
Berkeley). 



teresting proposals, e.g. for tunable spin filters are based 
on that fact f58ll . However, the existence of the edge states 
has been challenged by DFT calculations [59 1 revealing 
a thermodynamic instability of the edge states in realistic 
environments. STS measurements have been performed 
on nanoribbons, which are produced by zipping calcinated 
nanotubes in an organic solution and subsequently using 
an ultrasonic bath ll60l . The nanoribbons are dispersed on 
a Au(lll) surface. Figure [6| a) shows an STM image of 
such a nanoribbon. The edge area forms a little bump of 
about 4-5 A in height indicating bond formation of the 
edge with the substrate. From the atomic resolution on 
the bumped region, one can identify the edge geometry as 
(8,1) meaning that the edge consists of repeated segments 
of 8 zig-zag unit cells and a perpendicular armchair unit 
cell. The spectroscopy at the edge shown in Fig. |6|b) re- 
veals a double peak around the Fermi level in line with 
expectations from two magnetic edge states at the two 
ribbon edges, which are coupled antiferromagnetically to 
each other. The double peak disappears in the direction 
perpendicular to the edge after 2-3 unit cells. Figure 6]c) 
shows the energy distance A of the peaks measured for 
nanoribbons of different widths. The value of A decays 
with width as expected from the antiferromagnetic cou- 
pling of the edges. The values for A are in quantitative 



correspondence with a Hubbard model using a hopping 
term t = 2.7 eV and an on-site energy U = 1.35 eV ll6D . 
This is good evidence that the edge state can survive even 
in solution, although a direct proof of its magnetic charac- 
ter is desirable. 



5.3 Landau levels A nice fingerprint for STS sensi- 
tivity are Landau levels, which appear in magnetic field B 
due to the orbital quantization of electrons. In 2D systems, 
the quantization is complete leading to intricate transport 
effects as the integer quantum Hall effect ll62l or the frac- 
tional quantum Hall effect |63 1. Both effects have been ob- 
served on graphene IT1 I261 . The fractional quantum Hall 
effect, however, requires a low disorder potential, which 
has only been achieved, e.g., by freely suspending the 
graphene [26] or depositing it onto an intermediate layer 
of BN(0001) between a SiC>2 substrate and graphene |27|. 
If the disorder is reduced sufficiently, it is found that the 
thermal stability of the fractional plateaus is much larger 
than for the usual GaAs 2D systems [ 26 . 27 1 . This is ex- 
pected, since the electrons are vertically quenched into a 
single layer, while GaAs 2D systems exhibit thicknesses 
of several nm. Thus, the electron-electron interaction at the 
same electron density is significantly larger in graphene. 
Landau levels detected by STS, which have previously 
been observed on semiconductor 2D systems [64], have 
been first found for graphene layers on highly oriented py- 
rolytic graphite (HOPG) as shown in Fig. [TJa) [51 1. The 
graphene layer has a slightly larger distance to the under- 
lying graphene layer than within HOPG, which appears to 
be sufficient to electronically decouple the graphene layer 
from the substrate. The fingerprint of the graphene Landau 
levels is its energy dependence = ±ud • yj2efi ■ \N\B 
(«d — 10 6 m/s: Dirac velocity, e: electron charge, N € Z: 
Landau level index) [3[. The peaks shown in Fig.|5ja) in- 
deed exhibit such a behavior |5ll . Graphene Landau levels 
have also been observed on SiC(0001) [65 1, which is the 
carbon terminated face. The sublimation method leads to 
thicker graphite layers, however, with individual graphene 
layers rotated with respect to each other l66l . It turned out 
that the rotation leads to an effective decoupling of the in- 
dividual layers, resulting, e.g., in graphene cyclotron ab- 
sorbtion already at B = 0.2 T [67|. The corresponding 
Landau levels are shown in Fig. |7Jb) with the constituent 
^J\N\B dependence displayed in Fig.^c) 1651 . The Lan- 
dau levels have been used to map the disorder potential, 
which fluctuates on a length scale of about 50 nm by ±5 
meV [65 1 being very similar to the potential fluctuations 
found for graphene on BN (see Fig. [5]T)) or on InAs 2D 
systems 11681 . First attempts to measure Landau levels on 
graphene flakes on Si02 were not successful |50| because 
of the strong disorder potential which could fluctuate by 
±150 meV (see Fig. |5Je)). However, by chlorinating the 
substrate, the disorder is sufficiently reduced such that Lan- 
dau levels are observed (see Fig.|7Jl) l69l . By applying the 
gate voltage as shown in Fig.[7]e), one can nicely observe 
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Figure 7 (a) dl/dV curves recorded on graphene on HOPG at different B fields as marked, Landau level indices N are 
marked at the peaks of the (10 T)-curve, V stab = 0.3 V, igtab = 0.02 nA, T = 4.4 K EQ; (b) dl/dV curves recorded 
on graphene on top of several graphene layers on SiC(0001) at different B fields as marked , V B tah = 0.35 V, / sta b = 0.4 
nA, T = 4.3 K 1651 : (c) observed peak energies from (b) as a function of sgn(iV) y/\N\B (N € Z) ; (d) dl/dV curves 
recorded on graphene on chlorinated SiC>2 at different B fields as marked, V^tab = 0.3 V, 7 s tab = 0.02 nA, T = 4.4 K; 
(e) color plot of dl/dV intensity as a function of sample bias V and gate voltage Vg atc at B = 12 T; other parameters as 
in (d), Landau level index N is marked [69 1 ((a), (d), (e): courtesy of E. Andrei, Rutgers University, (b), (c): courtesy of J. 
Stroscio, NIST Gaithersburg). 



the jumps of Landau levels to the Fermi level, e.g. between 
N = -1 and N = around Vg ato = 30 V. 

5.4 Drift States The states corresponding to the Lan- 
dau levels are called drift states. They meander along 
equipotential lines of the disorder potential and exhibit a 
width of the magnetic length lg = y/h/ (eB) [25|. Most 
drift states are localized in potential valleys or at potential 
hills except of one state located in the energetic center of 
the Landau level. This state traverses the whole potential 
landscape and, thus, is called the extended state. It is re- 
sponsible for the quantum Hall transitions. The transition 
from localized to extended states and back to localized 
states has been imaged previously on InAs 2D systems 
ll70l . It has also been found for the graphene layer on 
SiC(000T) ED. Fig.[|:a)-(d) show dl/dV images at B = 8 
T taken at energies around the Landau level with index 
N — 0. The magnetic length is marked in Fig. [SJd) and 
indeed corresponds to the widths of the observed mean- 



dering structures. A transition from localized states within 
the potential valleys ((a),(b)) via an extended state ((c)) to 
localized states around the potential hills ((d)) is apparent. 
Experiments on graphene flakes on SiC>2 did not show 
drift states. Instead, the strong disorder leads to strongly 
localized charge puddles, which exhibit Coulomb block- 
ade peaks in STS experiments [72 1 . Chlorinating the SiC>2 
improves the disorder and indications of drift states be- 
come visible f73l . 

A more detailed analysis of the drift states on SiC(0001) 
reveals that a portion of the surface potential is given by the 
lattice mismatch and rotation angle of the probed graphene 
layer with respect to the underlying graphene layer. This 
implies a moire pattern and a periodic potential with a lat- 
tice constant of about 70 nm 1711 . Interestingly, the dl/dV 
curves exhibit a sublattice splitting in certain areas of this 
potential. Figure |8je) shows a line scan of dl/dV curves. 
In the central area, the peak energy jumps between two 
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Figure 8 (a)— (d) dl/dV images recorded on graphene on SiC(000T) at B — 8 T at the voltages marked, / = 0.4 nA, 
T = 4.3 K, magnetic length Is = 9 nm is marked in (d), I and II in (c) mark areas similar to the ones, where the histograms 
in (f) have been obtained; (e) color plot of dl/dV curves along a line showing the development of the N = Landau 
level, Vstab = 0.35 V, I s tab = 0.4 nA; I, II as in (c); (f) histogram of peak values of dl/dV curves as displayed in (e) fTTl ; 
(g) color plots of dl/dV curves around the Fermi level along a line displayed at different B fields as marked, spin levels 
and (K,K') levels are indicated at the top part, v is the local filling factor obtained from counting the number of peaks 
below the Fermi level (hole Landau levels are not observed for an unknown reason), half filling factors are interpolated 
and marked, since contributing a third level, T4tab = 0.25 V, I s tab = 0.2 nA, Vmod = 50 fi"V, T = 0.013 K [74] (courtesy 
of J. Stroscio, NIST Gaithersburg). 



values. Histograms of the peak values obtained in such 
areas are shown in the upper part of Fig. [8|f). They reveal 
two dominating peak values being separated by about 7 
meV. Other areas within the moire structure do not exhibit 
a peak splitting as demonstrated by the histogram in the 
lower part of Fig. |SJg). A plot of the peak splitting value 
as a function of position reproduces the moire pattern. The 
splitting is, thus, caused by the different stacking of the 
two graphene sublattices onto the underlying layer within 
these areas. Equivalently, one can say that the stacking of 
graphene to the sublayer locally lifts the valley degeneracy 
of the zero-order Landau level. 

At very low temperature (T = 13 mK), this system shows 
also spin splitting as demonstrated by the color plots in 
Fig. |8jg) |74|. For example, at B = 12 T, pairs of lines 
separated by the Zeeman energy of 1 .4 meV are observed 
above and below the Fermi level. Each pair of lines cor- 
responds to one of the sublattices at the measurement 



position. The filing factor v = 4 indicates that two more 
lines are observed below the Fermi level corresponding 
to the electron states of the N = Landau level. Thus, 
the four lines correspond to the N = 1 Landau level. 
By decreasing the magnetic field, the degeneracy of each 
line decreases according to till = eB/h («ll: number 
of states per m 2 in one spin- and valley-resolved Landau 
level). Thus, the levels above the Fermi level Ep must 
cross Ep, if the electron density does not change. This 
crossing is visible in Fig. [8jg). At the crossing point of 
the first spin level, the distance between the spin levels 
dramatically increases and it decreases again, if both spin 
levels have crossed Ep. The effect is well known as ex- 
change enhancement based on the fact that electrons with 
the same spin feel a weaker Coulomb repulsion than the 
ones with opposite spins [75 1. Moreover, a weak third line 
is observed at half-valued filling factor (marked by arrows 
in Fig. |HJg)), which is interpreted as an interaction effect 
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with the quantized 2D system of the underlying graphene 
layer. Notice that the exchange enhancement is the first 
tunable electron-electron interaction effect observed by 
STM on graphene. 

Experiments of graphene flakes on BN/SiC>2 are under way 
and one should be curious if signatures of the fractional 
quantum Hall plateaus [27 1 will be observed in STS. 

6 Conclusion Although a number of very interesting 
STM and STS results on graphene have already been ob- 
tained , there are still multiple questions to be solved. One 
of the most interesting ones is the appearance of fractional 
quantum Hall states, where very little is known about. 
Other interesting aspects include the magnetism of the 
graphene edge state [57 1, signatures of the so-called Zitter- 
bewegung [76] or the proposed charge instability at highly 
charged defects [23 1. While most results so far took advan- 
tage of the fact that a 2D system is directly at the surface, 
future experiments might explore the specifics of the Dirac 
nature of the particles. This might substantiate the idea to 
mimic elementary particle physics by small scale exper- 
iments and promises still an exciting time with STM on 
graphene. Most likely, also the various applications now 
arising for graphene might pose questions to be solved on 
the nano-scale. 
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